Impact of CLAS and COMPASS data on Polarized Parton Densities and Higher Twist 
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We have re-analyzed the world data on inclusive polarized DIS including the very precise CLAS 
proton and deuteron data, as well as the latest COMPASS data on the asymmetry Af, and have 
studied the impact of these data on polarized parton densities and higher twist effects. We demon- 
strate that the low CLAS data improve essentially our knowledge of higher twist corrections 
to the spin structure function gi, while the large COMPASS data influence mainly the strange 
quark density. In our new analysis we find that a negative polarized gluon density, or one that 
changes sign as a function of x, cannot be ruled out on the basis of the present DIS data. 
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I. INTRODUCTION 

The European Muon Collaboration (EMC) experiment 
[l[ at CERN found that a surprisingly small fraction of 
the proton spin is carried by the spin of the quarks. This 
observation was a big challenge to our understanding of 
the partonic spin structure of the nucleon, i.e., how the 
nucleon spin is built up out from the intrinsic spin and or- 
bital angular momentum of its constituents, quarks and 
gluons. Since that time substantial efforts, both exper- 
imental and theoretical, have been made to answer this 
question. Our present knowledge about the spin struc- 
ture of the nucleon comes mainly from polarized inclusive 
and semi-inclusive DIS experiments at SLAC, CERN, 
DESY and JLab, polarized proton-proton collisions at 
RHIC and polarized photoproduction experiments. One 
of the important and best studied aspects of this knowl- 
edge is the determination of the longitudinal polarized 
parton densities in QCD and their first moments 0, 
which correspond to the spins carried by the quarks and 
gluons in the nucleon. 

One of the features of polarized DIS is that a lot of 
the present data are in the preasymptotic region {Q^ ~ 
1 - 5 GeV^, 4 GeV2 < < 10 GeV^). This is es- 
pecially the case for the experiments performed at the 
Jefferson Laboratory. As was shown in to confront 
correctly the QCD predictions to the experimental data 
including the preasymptotic region, the non-perturhative 
higher twist (powers in l/Q^) corrections to the nucleon 
spin structure functions have to be taken into account 
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too. 

In this paper we study the impact of the recent very 
precise CLAS [5] and COMPASS [S] inclusive polarized 
DIS data on the determination of both the longitudinal 
polarized parton densities (PPD) in the nucleon and the 
higher twist (HT) effects. These experiments give im- 
portant information about the nucleon structure in quite 
different kinematic regions. While the CLAS data en- 
tirely belong to the preasymptotic region and as one can 
expect they should mainly influence the higher twist ef- 
fects, the COMPASS data on the spin asymmetry Af 
are large data and they should affect mainly the po- 
larized parton densities. In addition, due to COMPASS 
measurements we have for the first time accurate data 
at small x (0.004 < a; < 0.015), where the behaviour of 
the spin structure function gf should be more sensitive 
to the sign of the gluon polarization. 



II. NEXT TO LEADING QCD ANALYSIS OF 
THE DATA 



In QCD the spin structure fimction gi has the following 
form (Q^ >> A^): 

gi{x,Q'^) = gi{x,Q'^)i,T + gi{x,Q'^)BT , (1) 

where "LT" denotes the leading twist (t = 2) contri- 
bution to gi, while "HT" denotes the contribution to 
gi arising from QCD operators of higher twist, namely 
T > 3. In Eq. ([T]) (the nucleon target label N is dropped) 

gi{x,Q\T = gi{x,Q\QCB + h™''{x,Q^)/Q^ 

+0(mVQ') , (2) 
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where 51(2;, Q^)pQCD is the well known (logarithmic in 
Q2) nLO pQCD contribution 

and h'^^'~^{x,Q'^) are the calculable kinematic target 
mass corrections 7], which effectively belong to the LT 
term. In Eq. ([3])", Aq{x,Q^), Aqix,Q^) and AG{x,Q^) 
are quark, anti-quark and gluon polarized densities in 
the proton, which evolve in according to the spin- 
dependent NLO DGLAP equations. dC{x)q^G are the 
NLO spin-dependent Wilson coefficient functions and the 
symbol ® denotes the usual convolution in Bjorken x 
space. Nf is the number of active flavors (Nf = 3 in 
our analysis). In addition to the LT contribution, the 
dynamical higher twist effects 

g2)HT = h{x, Q^)/Q^ + o(AVg*) , (4) 

must be taken into account at low Q^. The latter are non- 
perturbative effects and cannot be calculated in a model 
independent way. That is why we prefer to extract them 
directly from the experimental data. The method used 
to extract simultaneously the polarized parton densities 
and higher twist corrections to 171 is described in Q . Ac- 
cording to this method, the gi/Fi and gi/Fi) data 
have been fitted using the experimental data for the un- 
polarized structure function Fi{x,Q^) 



a simple parametrization 

xAut, (x, Qq) — TjiiAuX " xUy[x, Qq) , 

xAdy{x,Ql) = ridAdx'^''xdy{x,Ql), 

xAs{x,QI) = jjs AsX^'^xsI^XjQI), 

xAG{x,Ql) = r^gAgx''^xG{x,Ql), (6) 

where on the RHS of ^ we have used the MRST99 (cen- 
tral gluon) 10] parametrizations for the NLO (MS) unpo- 
larized densities. The normalization factors Ai in ([6]) are 
fixed such that rji are the first moments of the polar- 
ized densities. The first moments of the valence quark 
densities rju and rjd are constrained by the baryon decay 
constants (F-fD) and (3F-D) Q assuming SU(3)f sym- 
metry. Bearing in mind that the light quark sea densi- 
ties Au and Ad cannot, in principle, be determined from 
the present inclusive data (in the absence of polararized 
charged current neutrino experiments) we have adopted 
the convention of a flavor symmetric sea 

AUsea = Aw = Adsea ^ Ad ^ As = As. (7) 

Note that this convention only affects the results for 
the valence parton densities, but not the results for the 
strange sea quark and gluon densities. 

In polarized DIS the range and the accuracy of the 
data are much smaller than that in the unpolarized case. 
That is why, in all calculations we have used a fixed value 
of the QCD parameter Aj^(n/ = 4) = 300 MeV, which 
corresponds to as(M|) = 0.1175, as obtained by the 
MRST NLO QCD analysis [IH of the world unpolarized 
data. This is in excellent agreement with the current 
world average as(M2) = 0.1176 ± 0.002 pj]. 
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(5) 
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As usual, Fi is replaced by its expression in terms of 
the usually extracted from unpolarized DIS experiments 
F2 and R and phenomenologicalparametrizations of the 
experimental data for F2{x, Q^) [8] and the ratio R{x, Q^) 
of the longitudinal to transverse 7iV cross-sections [9| are 
used. Note that such a procedure is equivalent to a fit 
to (51) exp, but it is more precise than the fit to the gi 
data themselves actually presented by the experimental 
groups because here the gi data are extracted in the same 
way for all of the data sets. Note also, that in our analysis 
the logarithmic dependence of h{x,Q-^) in Eq. ([5]), 
which is not known in QCD, is neglected. Compared to 
the principal l/Q^ dependence it is expected to be small 
and the accuracy of the present data does not allow its 
determination. Therefore, the extracted from the data 
values of h{x) correspond to the mean for each x- 
bean (see Table II and the discussion below). 

As in our previous analyses, for the input NLO polar- 
ized parton densities at Qq = 1 GeV^ we have adopted 



III. RESULTS OF ANALYSIS 

In this section we will discuss how inclusion of the 
CLAS proton and deuteron gi/Fi data 0] and the new 
COMPASS data on Af Q influence our previous results 
[1] on polarized PD and higher twist obtained from the 
NLO QCD fit to the world data P, [H, [11, before the 
CLAS and the latest COMPASS data were available. 



A. Impact of CLAS data 

The CLAS EGi/p, d data (633 experimental points) 
we have used in our analysis are high-precision data in 
the following kinematic region: {x ^ 0.1 — 0.6, ~ 
1 - 5 GeV^, W > 2 GeV}. As the CLAS data are 
mainly low data where the role of HT becomes im- 
portant, they should help to fix better the higher twist 
effects. Indeed, due to the CLAS data, the determina- 
tion of HT corrections to the proton and neutron spin 
structure functions, hP{x) and h"(x), is significantly im- 
proved in the CLAS x region, compared to the values 
of HT obtained from our LSS'05 analysis 3 in which a 
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FIG. 1: Effect of CLAS data on tlie liiglier twist values. 



NLO(MS) QCD approximation for gi{x, Q^)l,t was used 
(see Table I). This effect is illustrated in Fig. 1. One 
can conclude now that the HT corrections for the proton 
target are definitely different from zero and negative in 
the X region: 0.1-0.4. Also, including the CLAS data in 
the analysis, the HT corrections for the neutron target 
are better determined in the x region: 0.2-0.4. Note that 
h"'{x) at a; 0.5 was already fixed very precisely from 
the JLab Hall A data on the ratio g'^'' /f["'\ 

The values obtained for the parameters of the input po- 
larized PD are presented in Table I and compared with 
those of LSS'05. Note that the extracted polarized PD 
correspond to the Set 2 of NLO(MS) LSS'05 PPD. As 
expected, the central values of the polarized PD are prac- 
tically not affected by the CLAS data (see Table I). This 
is a consequence of the fact that at low the devia- 
tion from logarithmic in pQCD behaviour of gi is 
accounted for by the higher twist term ^ in gi. Indeed, 
if one calculates the x^-probability for the combined 
world-l-CLAS data set using the LSS'05 polarized PD 
and corresponding HT values, the result for is 938.9 
for 823 experimental points, which significantly decreases 
to 718.0 after the fit. As seen from Table I, the best 
fit to the combined data is achieved mainly through the 
changes in the HT values. This supports the theoretical 
framework in which the leading twist QCD contribution 
is supplemented by higher twist terms of ©(Aqj-^q/Q^). 
One can see also from Table I, that the accuracy of the 
determination of polarized PD is essentially improved. 
This improvement (illustrated in Fig. 2) is a consequence 
of the much better determination of higher twist contri- 
butions to the spin structure function gi, as discussed 
above. Due to the good accuracy of the CLAS data, one 
can split the measured x region of the world-t-CLAS data 
set into 7 bins instead of 5, as used up to now, and there- 



TABLE I: The parameters of the NLO(MS) input PPD at 
= 1 GeV^ and HT as obtained from the best fits to the 
world data (LSS'05) and combined world-hCLAsH 

data set (LSS'06). The errors shown are total (statistical and 
systematic). The parameters marked by (*) are fixed. 



Fit 


LSS'05 (Set 2) 


LSS'06 


DF 


190 - 16 


823 - 16 




154.5 


718.0 


X'/DF 


0.888 


0.890 




0.926* 


0.926* 




0.252 ± 0.037 


0.252 ± 0.025 


Vd 


- 0.341* 


-0.341* 




0.166 ± 0.124 


0.166 ± 0.092 


■ns 


- 0.070 ± 0.008 


- 0.070 ± 0.007 


as 


0.656 ± 0.069 


0.679 ± 0.046 


Vg 


0.179 ± 0.267 


0.296 ± 0.197 


ag 


2.218 ± 1.650 


2.465 ± 0.878 





h^ix,) 


[GeV 


'] 


0.028 


0.018 ± 0.047 




0.004 ± 0.040 


0.100 


- 0.031 ± 0.032 




0.049 ± 0.013 


0.200 


- 0.100 ± 0.040 




0.045 ± 0.012 


0.350 


0.004 ± 0.046 




0.041 ± 0.011 


0.600 


0.036 ± 0.020 




0.023 ± 0.013 




h"{x,) 


[GeV 




0.028 


0.182 ± 0.065 




0.168 ± 0.056 


0.100 


0.196 ± 0.038 




0.197 ± 0.033 


0.200 


0.081 ± 0.061 




0.047 ± 0.029 


0.325 


0.025 ± 0.029 




0.025 ± 0.020 


0.500 


0.014 ± 0.013 




0.015 ± 0.011 



fore, can determine more precisely the x-dependence of 
the HT corrections to gi. The numerical results of the 
best fit to the data using 7 x-bins are listed in Table II 
(first column). In Fig. 3 the HT values corresponding 
to 5 and 7 x-bins are compared. As seen in Fig. 3, the 
more detailed x-space behaviour of the HT contribution, 
obtained when using 7 x-bins, suggests a smoother func- 
tion dependence in x and will help us to calculate more 
precisely their first moments in the experimental x re- 
gion and to compare them with the predictions given by 
different models. This point will be discussed below. 

Comparing the fitted PPD parameters corresponding 
to 5 and 7 x-bins for HT (Table I, 2nd column and Table 
II, 1st column, respectively), we observe that the input 
valence Au^ and Ad^, as well as the strange quark sea As 
densities are practically identical. The only exception is 
the gluon density, but as seen in Fig. 4, the curve of the 
gluon density corresponding to 7 x-bins lies within the 
error band of a;AG(5 bins). The curves corresponding 
to (Au -I- Au) and (Ac? -I- Ad) densities cannot be distin- 
guished in the experimental region and for that reason 
they are not shown in Fig. 4. Note that the curves 
corresponding to input As{x) at Qo = 1 GeV^ also can- 
not be distinguished. They became slightly different at 
Qo (^^^ Pi§- 4) because of the mixture between 
the gluons and sea quarks due to their evolution. We 
consider that the small correlation between the polarized 
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FIG. 2: Impact of CLAS data on the uncertainties for NLO(MS) polarized quark and gluon densities. 
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FIG. 3: Comparison between the higher twist values corre- 
sponding to 5 and 7 a;-bins. 



gluon density and the HT corrections we have found re- 
flects the fact that the gluons are not well constrained 
from the present inclusive DIS data. 
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FIG. 4: Comparison between polarized NLO(MS) LSS'06 
strange quark and gluon densities corresponding to fits of the 
data using 5 and 7 x-bins for higher twist. 



B. Impact of new COMPASS data 



When this analysis was finished, the COMPASS Col- 
laboration at CERN reported new data on the longitu- 
dinal asymmetry A\ Q. The new data are based on a 
2.5 times larger statistics than those of Ref. [3| used in 
our analysis. In contrast to the CLAS data, the COM- 
PASS data are at large and are the only precise data 



covering the low x region: 0.004 < x < 0.015, where 
the behaviour of the spin structure function g'^ should be 
more sensitive to the sign of the gluon polarization. Note 
also, that due to the larger statistics the latest COM- 
PASS data give more precise and detailed information 
about A\ and in the above experimental region (see 



5 





-0.2 

(b) 



new COMPASS 
■ ■ (AG > 0) 
-- g/T (AG>0) 
■■ g,^^ (AG<0) 



0.1 



-0.025 





— Z.J vjc V 




COMPASS 






- new data 






old data 




xAs ■ 



Q' = 2.5GeV' 



xAG 




FIG. 6: Effect of new COMPASS data on the NLO(MS) 
LSS'06 polarized parton densities. 



FIG. 5: Comparison of our NLO(MS) results for Ai (a) and 
gf (b) corresponding to AG > and AG < with the new 
COMPASS data at measured x and values. Error bars 
represent the total (statistical and systematic) errors. 



Fig. 5). 

In view of this, we have re-analyzed the data super- 
seding the old set of COMPASS data with the latest one 
in order to study the impact of the new COMPASS data 
on the results reported above. 7 a;-bins for extracting 
the HT values were used in the fits. The numerical re- 
sults are listed in Table II. As mentioned in Section II, 
the logarithmic dependence of higher twist Q^) is 
neglected in our analysis. So, the numerical values h{xi) 
presented in Table II correspond to the mean value of Q'^ 
for any .T-bin. Keeping in mind that the higher twist val- 
ues h{xi) arc mainly determined from the preasymptotic 
region (Q^ _ 1 - 5 GcV^, W > 2 GeV), the mean val- 
ues < Q? > in Table II correspond to the experimental 
points for each x-bin with <h GeV^. 

The QCD theoretical curves for Af corresponding to 
the best fits to the data with positive and negative gluon 
polarization are shown in Fig. 5(a). The old COMPASS 
data and Af calculated using the LSS'06 polarized PD 
and HT corrections (discussed in the previous Section), 
are also presented. Note that for a; > 0.1 the theoretical 
curves corresponding to the fits using the new or old set 
of the COMPASS data cannot be distinguished and for 



that reason this x region is not shown in Fig. 5(a). The 
best fit to the new gi data is illustrated in Fig. 5(b). 

The effect of the new data on the polarized parton 
densities and the higher twist corrections is illustrated 
in Fig. 6 and Fig. 7, respectively. While (Au + Au) 
and {Ad + Ad) parton densities do not change in the 
experimental region (for that reason they are not shown 
in Fig. 6), the magnitudes of both the polarized gluon 
and strange quark sea densities and their first moments 
slightly decrease (see Fig. 6 and Table II). As a con- 
sequence, AS((5^ = 1 GeV^) increases from (0.165 ± 
0.044) to (0.207 ± 0.040) for AG > and (0.243 ± 0.065) 
for AG < (sec below the discussion about AG < 0). 

As the COMPASS data are mainly at large Q^, the im- 
pact of the new data on the values of higher twist correc- 
tions is negligible, and as expected, they do not improve 
the uncertainties of HT. The new central values practi- 
cally coincide with the old ones (see Fig. 7(a)). The only 
exception are the central values of HT at small x for both 
the proton and the neutron targets which are slightly 
lower than the old ones. Note that this is the only region 
where the COMPASS DIS events are at small Q^. 1.4 
GeV^. As a result, in the small x region two opposite ten- 
dencies occur. In order to make gf consistent with zero 
for X < 0.03, the HT contribution /i'* = {hP + h"-)Q. 925/2, 
which is positive, decreases slightly, while (s'i)lt, which 
is negative, grows slightly due to the smaller negative 
contribution of As{x, Q^) and the smaller contribution of 
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TABLE 11: Effect of tlie new COMPASS data on polarized PD and HT. Tlie parameters of tfie NLO(MS) input PPD are given 
at = 1 GeV'^. 7 bins in X liave been used to extract tlie HT values. The errors shown are total (statistical and systematic). 



The parameters 


marked by (*) 


are fixed. 






Fit 




LSS'06 [old COMPASS) 


LSS'06 [new 


COMPASS) 






AG > 


AG > 


AG < 


DF 




823 - 20 


826 - 20 


826 - 20 






716.2 


721.7 


722.9 






0.892 


0.895 


0.897 






0.926* 


0.926* 


0.926* 


au 




0.264 ± 0.027 


0.273 ± 0.028 


0.273 ± 0.028 


Vd 




- 0.341* 


-0.341* 


-0.341* 


Cld 




0.172 ± 0.118 


0.202 ± 0.118 


0.160 ± 0.108 


Vs 




- 0.070 ± 0.006 


- 0.063 ± 0.005 


- 0.057 ± 0.010 


ds 




0.674 ± 0.053 


0.715 ± 0.052 


0.746 ± 0.088 


Vg 




0.173 ± 0.184 


0.129 ± 0.166 


- 0.200 ± 0.414 


ttg 




2.969 ± 1.437 


3.265 ± 1.668 


0.698 ± 0.806 


Xi 


<Qt> 




h^ix^) [GeV] 




0.028 


2.0 


0.034 ± 0.040 


0.010 ± 0.039 


0.017 ± 0.041 


0.075 


2.4 


- 0.001 ± 0.030 


- 0.016 ± 0.030 


- 0.019 ± 0.037 


0.150 


1.7 


- 0.046 ± 0.010 


- 0.050 ± 0.009 


- 0.056 ± 0.018 


0.250 


1.8 


- 0.055 ± 0.011 


- 0.059 ± 0.010 


- 0.067 ± 0.013 


0.350 


2.4 


- 0.050 ± 0.013 


- 0.054 ± 0.012 


- 0.060 ± 0.013 


0.450 


3.2 


- 0.012 ± 0.015 


- 0.016 ± 0.015 


- 0.020 ± 0.015 


0.625 


4.1 


0.018 ± 0.015 


0.016 ± 0.015 


0.014 ± 0.015 


Xi 


< > 




h"{xi) [GeV^] 




0.028 


1.8 


0.187 ± 0.064 


0.165 ± 0.064 


0.180 ± 0.065 


0.075 


2.4 


0.190 ± 0.044 


0.173 ± 0.044 


0.174 ± 0.049 


0.150 


1.4 


0.104 ± 0.041 


0.107 ± 0.039 


0.092 ± 0.040 


0.250 


1.5 


0.018 ± 0.030 


0.019 ± 0.030 


0.006 ± 0.029 


0.350 


2.2 


0.028 ± 0.026 


0.031 ± 0.025 


0.019 ± 0.023 


0.450 


3.0 


0.010 ± 0.020 


0.013 ± 0.020 


0.005 ± 0.019 


0.625 


3.9 


0.014 ± 0.013 


0.016 ± 0.013 


0.012 ± 0.012 



AG{x,Q^), convoluted with its Wilson coefficient func- 
tion 5Cg{x), which is negative in this x range (see Eq. 
®)- 

We have also checked the stability of our results with 
respect to a change in as(M^), which in our analysis 
coincides with its current world average, as mentioned 
in Section II. When as{M^) is varied by one standard 
deviation ±0.002, the change of the values of the free pa- 
rameters is within their errors. In particular, the change 
of r/s and rjg, the first moments of the polarized quark sea 
and gluon densities, is smaller than 10% of their standard 
deviations. 



C. The sign of the gluon polarization 

We have also studied the possibility of a negative po- 
larized gluon density. Starting with a negative value for 
Tjg — AG'(Qq) (the first moment for the input gluon po- 
larized density AG(x,(3§)), we have found a minimum 
in corresponding to a negative solution for r]g, and 
to negative AG^Q"^) and xAG^x^Q"^). The values of 
corresponding to the fits with AG > and AG < are 
practically the same (see Table II) and the data cannot 
distinguish between these two solutions for AG (see Fig. 



5(a)). Note that in our previous analyses we also found 
solutions with negative AG, but they were not presented 
because the corresponding were significantly larger 
than those corresponding to the solutions with positive 
AG. 

In Fig. 8 the negative polarized gluon density is com- 
pared with the positive one. As seen, the shape of the 
negative gluon density differs from that of positive one, 
but in both cases the magnitude of a;AG is small. Con- 
sequently the parton densities obtained in the fits with 
AG > and AG < are almost identical. For the 
strange quarks this is illustrated in Fig. 8. Thus the the- 
oretical curves (5i )lt for the two types of gluon polariza- 
tion are practically identical, even in the region x < 0.01 
(see Fig. 5(b)). 

Furthermore, as seen in Fig. 7(b), the extracted HT 
values corresponding to the fits with AG > and AG < 
0, are effectively identical. Thus also the total theoretical 
expression {gf)tot is essentially the same for AG > and 
AG < 0, even at very small x < 0.01. 

These results are in contrast to those obtained in the 
COMPASS analysis Q where there is a significant differ- 
ence between the theoretical curves corresponding to the 
cases AG > and AG < at very small x, i.e. in the 
region 0.004 < x < 0.02. The reason for this lies in the 
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FIG. 7: Effect of new COMPASS data on the higher twist values (a). Comparison between HT values corresponding to the fits 
with AG > and AG < (b). 
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FIG. 8: Strange quark sea densities xAs(x) corresponding to 
the fits with AG > 0, AG < and changing in sign xAG. 



question of HT contributions, which are not taken into 
account by COMPASS. In the above x region, is small 
(Q^ ^ 1 _ 3 GeV^) and we have found that the HT con- 
tribution to (3i)tot, h'^{x)/Q^, is positive and large, up 
to 40% of the magnitude of (gf)LT (see Fig. 5(b)). Thus 
what is fitted by (g5*)LT (COMPASS) is significantly dif- 
ferent from what is fitted by our ((7f)LT(LSS) at small x, 



i.e. (3^)lt(C0MPASS) = (gf)LT(LSS)+hd(x)/Q2. As 
a result: i) The strange quark sea densities obtained in 
the two analyses are different, especially in the case of 
AG < (see Fig. 9 and Fig. 10). ii) The gluon den- 
sities obtained by COMPASS in both fits (AG > and 
AG < 0) are more peaked than ours. 

Finally, concerning the possible solution with negative 
AG we would like to point out the much larger uncertain- 
ties in the determination of the strange quark sea and 
gluon densities, xAs and xAG, and respectively, their 
first moments (see Fig. 11 and Table II). As seen from 
Fig. 11, the positive gluon density xAG{x) lies in the er- 
ror band of the negative gluons except for x larger than 
0.2. xAs{x) corresponding to the positive AG solution 
lies entirely in the error band of a;As(AG < 0). 

Bearing in mind the high precision of the CLAS and 
new COMPASS data over a large range in we have 
studied the possibility to obtain from the fit to the world 
inclusive DIS data a gluon density which changes sign 
as a function of x. Such a density was discussed in [1^ 
in order to describe the double longitudinal spin asym- 
metry All of inclusive tt" production in polarized P+P 
collisions measured by the PHENIX pj] and STAR ^ 
Collaborations at RHIC. To that end we introduced a fac- 
tor (1 -1-70;*) in the input gluon density in ^ with two 
new free parameters, 7 and S, to be determined from the 
fit to the data. In Fig. 8, the determined strange quark 
and gluon densities at — 2.5 GeV^ are compared with 
those corresponding to the positive and negative AG so- 
lutions. As seen from Fig. 8, the oscillating in sign gluon 
density lies between those of positive and negative AG. 
The value of per degree of freedom is 0.895, which 
coincides with the values obtained with purely positive 
or negative a;AG(a;). 
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FIG. 9: Comparison between our strange quark sea and gluon 
densities corresponding to AG > and those obtained by 
COMPASS [Ol. 



FIG. 11: The uncertainties for the strange quark sea and 
gluon densities corresponding to a negative gluon polariza- 
tion. 
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FIG. 10: Comparison between our strange quark sea and 
gluon densities corresponding to AG < and those obtained 
by COMPASS 0. 



Thus, we are forced to conclude that the accuracy and 
range of the present DIS data is not good enough 
to discriminate between these three possibilities. At 
= 1 GeV^, the shape of the oscillating in sign po- 
larized gluon density is consistent with that obtained by 
the AAC Collaboration from a combined analysis of DIS 
(CLAS and new COMPASS not included) and 7r° asym- 
metry data Note, however, that compared to the 
central value of the first moment AGaac — —0.56 ± 2.16 
at = 1 GeV"^, presented in fl8|, the central value of 
our AG is positive, 0.006, and much smaller in magni- 
tude. Under evolution in neither AG{Q'^)aac, nor 
our AG((5^) changes sign, and their magnitudes increase 
with increasing of Q^. As a result, the shape of the cor- 
responding gluon densities for > Qg will follow dif- 
ferent tendencies:a;AG(a;, (5^)aac becomes negative for 
larger x with increasing of Q^, while our gluon density 
for > 6 GeV^ is positive for any x in the experimental 
region (see Fig. 12). 

In Fig. 13 the ratio AG{x)/G{x) calculated for 
the different AG(a;) obtained in our analysis and using 
G(a;)MRST'02 taken from [l^, is compared to the exist- 
ing direct measurements of AG/G [20|. (Note that the 
MRST'02 unpolarized parton densities were used also in 
the positivity constraints imposed on the polarized par- 
ton densities obtained in our analysis.) The theoretical 
curves are given for = 3 GeV^. The most precise 
value for AG/G, the COMPASS one, is well consistent 
with any of the polarized gluon densities determined in 
our analysis. 
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Our values for the first moments for the proton, neu- 
tron and (Ji^ — h"') are consistent within the errors with 
those extracted directly from the analysis of the first mo- 
ments of and given in Rfs. [l^], [23'| and [24*1, respec- 
tively. Note that our value for the non-singlct (h'P — h^) 
is in agreement with the QCD sum rule estimates 12511 
as well as with the instanton model predictions [2^ |27| |. 
The values obtained for the non-singlet (ft^ — /i") and 
singlet (hP + h") quantities are in qualitative agreement 
with the relation \hP + h'^\ « \hP - h'^\ derived in the 
large hmit in QCD 



FIG. 12: Evolution in of oscillating-in-sign gluon density. 
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FIG. 13: Comparison between the experimental data and 
NLO(MS) curves for the gluon polarization AG(x)/G(a;) at 
= 3 GeV^ corresponding to AG > 0, AG < and an 
oscillating-in-sign lAG. Error bars represent the total (sta- 
tistical and systematic) errors. The horizontal bar on each 
point shows the a;-range of the measurement. 

Finally, let us briefly discuss the values of the first mo- 
ments of the higher twist corrections to the proton and 
neutron structure function g\. Using the values for {x) 
corresponding to AG > (second column in Table II) we 
obtain for their first moments in the experimental region: 

^0.75 

= \ h^{x)dx, {N=p,n) (8) 

Jo. 0045 

hP = (-0.014 ± 0.005) GeV^ for the proton and /i" = 
(0.037 ±0.008) GeV^ for the neutron target. As a resuh, 
for the non-singlet {hP — h") and the singlet {hP + h") we 
obtain (-0.051 ±0.009) GeV^ and (0.023 ± 0.009) GeV^, 
respectively. The errors are total (statistical and system- 
atic). The systematic errors are added qudratically. Note 

that in our notation h = h{x)dx = 4M^((i2 ± /2)/9, 
where 6,2 and /2 are the well known quantities, connected 
with the matrix elements of twist 3 and twist 4 operators, 
respectively 



IV. SUMMARY 

We have studied the impact of the CLAS and latest 
COMPASS data on the polarized parton densities and 
higher twist contributions. It was demonstrated that the 
inclusion of the low CLAS data in the NLO QCD anal- 
ysis of the world DIS data improves essentially our knowl- 
edge of higher twist corrections to the spin structure 
function gi. As a consequence, the uncertainties in the 
longitudinal polarized parton densities become smaller. 
The central values of the densities, however, are not af- 
fected and they practically coincide with those of LSS'05 
polarized parton densities determined from our previous 
analysis. In contrast to the CLAS data, the new more 
precise COMPASS data infiuence the strange quark den- 
sity, but practically do not change the HT corrections. 
Given that the COMPASS data is mainly at large Q^, 
this behaviour supports the QCD framework, in which 
the leading twist pQCD contribution is supplemented by 
higher twist terms of 0{K^(^j^/Q'^). 

We have observed that the fit to the world gi data 
involving the CLAS and new COMPASS data yields 
three possible solutions for the polarized gluon density, 
AG{x) > 0, AG{x) < and an changing-in-sign AG{x), 
which equally well describe the present DIS data. Also, 
all of them are in a good agreement with the directly 
measured quantity AG{x)/G{x) reported by COMPASS, 
although their shapes are very different. We have found 
that the magnitude of the gluon polarization is small, 
lAGI < 0.3 at Q2 ^ 1 GeV2. We have also found 
that the higher twist contribution to gf in the x range 
0.004 < X < 0.03 is positive and large, up to 40% of the 
magnitude of {gf)ur at {x = 0.0046, = 1.1 GeV^} 
and therefore, gf is not too sensitive to the sign of the 
gluon polarization in the above x region, when the higher 
twist corrections are taken into account. 
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